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Abstract 

A formalism based on the complex-scaling method is presented to solve the few particle scattering 
problem in configuration space using bound state techniques with trivial boundary conditions. 
Several applications to A=3,4 systems are presented to demonstrate the efficiency of the method 
in computing elastic as well as break-up reactions with Hamiltonians including both short and 
long-range interaction. 
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I. INTRODUCTION 



The theoretical description of the quantum-mechanical collisions is one of the most com- 
plex problems in theoretical physics. The main difficulty to solve the scattering problem in 
configuration space is related to the fact that, unlike the bound state case, the scattering 
wave functions are not localized. One is therefore obliged to solve multidimensional integro- 
differential equations with extremely complex boundary conditions. Therefore, finding a 
method which could enable us to solve the scattering problem without an explicit use of the 
asymptotic form of the wave function is of great importance. Recently, the research in this 
field has been intensified very much IrtZ]. 

One of the pioneering works in this direction is the development of the Lorentz integral 
transform l|, which allows to calculate the integral cross section of a scattering process using 
bound state like techniques. Still its application to compute differential observables becomes 
rather involved. Few years later, the complex energy method has been applied to the four- 
nucleon problem with separable potential [2], providing full information on the scattering 
process. Lately, after some technical improvements, this method was succesfully applied 
to describe realistic four-nucleon break-up process [7fl. Other recent developments include 
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momentum lattice technique [3] and a method based on discretized continuum solutions [6(, 
although they are no yet tested in the four-body sector and to the break-up case respectively. 

On the other hand, already in the late sixties, Nuttal and Cohen [8| proposed a very 
handy method to treat the scattering problem for short range potentials using coordinate 
space variational bound state techniques, namely the complex scaling method. However 
application of this method to the scattering problem has lingered. Only recently a vari- 
ant of the complex scaling method to calculate scattering observables above the break-up 
threshold has been proposed and applied by Giraud et al. [9|. Still this variant relies on the 
spectral function formalism and requires a diagonalization of the full N-body matrices to 
get converged results, which makes it difficult to extend beyond the N=3 case. 

In [ij] we have demonstrated that only a slightly modified original version of the complex 
scaling method by Nuttal and Cohen 8], if combined with Fadddeev equations, turns to be 
very efficient in describing three-nucleon scattering, including break-up reactions. We have 
shown that this method may treat strong interaction of almost any complexity: realistic 
local and non-local potentials, optical potentials, all of them in conjunction with repulsive 
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Coulomb force. The method allows to calculate both elastic as well as break-up amplitudes, 
thus providing a full description of the three-body collisions. There are no formal obstacles 
in extending this method to treat collisions involving any number of particles, as long as 
one is able to handle the eventual large scale numerical problem. In this contribution we 
will summarize this current effort and present some new results concerning the solution of 
the 4-nucleon scattering problem above the breakup thresholds. 



II. FORMALISM 



In our previous paper 5j we have presented the basic ideas of the complex scaling method 
used to solve the scattering problem for 2- and 3-particle systems. Here, without loosing 
generality, we will briefly summarize the formalism for solving the Faddeev-Yakubovski 
(FY) equations for a system of four identical particles. In this particular case, we start by 
separating the incoming plane wave from the FY components (see figured]): 

K(x, y, z) = K° u \x, y, z) + K in (x, y, z) (1) 
H(x, y, z) = H out (x } y, z) + H in (x, y, z). (2) 

Note, that the four-body FY equations involve components of two types, K and H, in the 
asymptote describing the elastic 3+1 and 2+2 particle channels respectively. Therefore, in 
the last equation, one has H m = by considering collision of 3+1 particle clusters whereas 
K m = by considering 2+2 particle collisions. Asymptotes of the FY components K out and 
H out will contain only various combinations of the outgoing waves. These components are 
solutions of the driven FY equations, which have the form: 

{E-H - V 12 ) K out - V 12 {P + + P~) [(1 + Q)K out + H out 

= V 12 (P+ + P~) \{l + Q)H m + QIC 
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It is interesting to introduce the complex scaling operator (CSO) 
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e dr 



8y ^ S 
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(3) 
(4) 

(5) 



where r 2 = x 2 + y 2 + z 2 and 9 is the complex scaling angle, of free choice. The action of the 
complex scaling operator on outgoing wave gives 



S exp(ikr) = exp(— kr sin 6) exp(ikr cos i 



(6) 





FIG. 1: The FY components Kf 2 3 and for a given particle ordering. As z — > oo, the K 
components describe 3+1 particle channels, while the H components contain asymptotic states of 
2+2 channels. 

Therefore if one acts on K out (or H out ) with the complex scaling operator (j5J, by keeping the 
complex scaling angle in the range < 9 < n/2, one gets exponentially decreasing functions 
K out = SK out and H out = SH out . By acting with the same operator on equations (EH), one 
gets equivalent equations for the scaled FY components K out and H outl : 



S(E-H - V l2 ) S~ l K° ut - SV l2 S-\P + + P~) [(1 + Q)K^ ut + lF ut 

= SV l2 S- 1 (P + + P-) [(1 + Q)W n + QK^ 
S(E-H - V l2 ) S- l FP ut - SV 12 S- l P [(1 + Q)K^ ut + P7^ ut 



SVioS^P 



1 + Q)K ir ' 



(7) 



(8) 



The functions to determine in the last equation, notably K out and H out , are exponentially 
bound. Thus, in principle, bound-state techniques based on square integrable basis func- 
tions, may be employed to solve these equations. Nevertheless, one should still ensure that 
the inhomogeneous terms on the right hand side of eqs. ([THE]) is also square integrable. 
Indeed, after the complex scaling, the incoming waves K in and H m , diverge when the inter- 
cluster separation increases. This divergence should be screened by the short range potential 
term SV^S -1 . Nevertheless, the direction in which propagates the incoming wave does not 
coincide completely with the direction of the potential term. For a system of four identical 



1 Note, that the four-body incoming wave components K m (or K m ) and if 1 " (or H ln ) are constructed 
from the corresponding bound-state wave functions of 3- and 2-particle systems respectively. 
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particles, the inhomogeneous terms are exponentially bound if the following conditions are 
satisfied: 

tan£<^, (9) 

in case of incoming wave of 1 + 3-type - B% being the binding energy of 3-particle cluster - 
and 

tan0<J|^, (10) 

V *- J crn 

when considering incoming wave of 2 + 2-type with a binding energy of 2-particle cluster 
equal B 2 . E cm denotes the energy of the two scattering clusters in the center of mass frame. 
Similar conditions for a 3-body problem have been derived in jsj. 

The scattering amplitudes are extracted from the solutions of eqs. (J3EJ) by using integral 
relations, obtained through the Green's theorem 

method described in our previous works (5), [ll|, 12j. Namely, first spin, isospin and angular 



lOj . Equations (J3EJ) are solved using the 



dependence of the FY components is expanded using partial-waves. The radial dependence 
of the FY components, in variables (x,y,z), is expanded by means of cubic spline basis and 
vanishing boundary conditions are implemented at the borders of the chosen 3-dimensional 
grid. 

III. RESULTS 

We will present in this section some chosen results obtained by applying the complex 
scaling method to 3- and 4-body systems. In table U we summarize the phase shifts and 
inelasticit y p arameters calculated for S-wave nucleon-deuteron scattering, using MT I-III 



reference 
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potential 13 1. Very nice agreement is obtained compared to the benchmark calculation of 
both by neglecting (n-d) as well as including (p-d) repulsive Coulomb 



interaction. Differential break-up amplitudes may be equally extracted with very nice accu- 
racy [5], although they are not reported in this work. 

More recently, calculations in 3-body sector have been extended to the study of reactions 
including compound nucleus. Here we present the results concerning d + 12 C — > p + 13 C 
scattering, which were obtained within a three-body model containing n, p and 12 C clusters. 
The interaction between the free nucleons is described by the realistic AVI 8 potential flJ |. 
whereas the iV — 12 C interaction is described by the optical CH89 potential [15:]. The n- 12 C 
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TABLE I: Phase shifts and inelasticity parameters calculated for S-wave nucleon-deuteron scatter- 
ing, using MT I-III potential [lj] 



n-d p-d 
P.W. E lah (MeV) This work Ref. [16, 17] This work Ref. [16] 



2 c 
b l 

2 


1/11 

14.1 


Ke(d) 


1U5.5 


1 n .in 

105.49 


i no a 

1U8.4 


1 no A i Tol 

108.41 [6\ 






V 


0.465 


0.4649 


0.498 


0.4983[1] 


2 


14.1 


Re{5) 


69.0 


68.95 


72.6 


72.60 






n 


0.978 


0.9782 


0.983 


0.9795[1] 


2 


42.0 


Re(5) 


41.5 


41.35 


43.8 


43.68[2] 






V 


0.502 


0.5022 


0.505 


0.5056 


2 


42.0 


Re{5) 


37.7 


37.71 


40.1 


39.96 [1] 






V 


0.903 


0.9033 


0.904 


0.9046 



potential in the 2 Pi partial wave is made real and is adjusted to support the ground state 
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of 13 C with 4.946 MeV of binding energy; the parameters are taken from Ref. 
adjustment is made in order to reproduce p — 13 C threshold in the transfer reactions. These 
results are summarized in Fig. [2] and compared to the ones obtained by the Lisboa group 19], 
who used the formalism of Alt-Grassberger-Sandhas (AGS) equations in momentum space. 
One may see an excellent agreement between these two calculations for d + 12 C and p + 13 C 
elastic cross sections, as well as for the transfer d + 12 C — > p + 13 C cross section. 

Concerning the four-body case, we present in table [II] a calculation of the integrated 
elastic and break-up cross sections for n — 3 H collisions. These results are obtained for the 
total isospin T = 1 channel, using MT I-III potential [13]. A rather reasonable agreement is 
obtained with the experimental values, although the agreement worsens at higher energies. 
MT I-III potential, being restricted to the S-waves, is not appropriate to this high energy 
domain. In figure [3] the corresponding differential elastic cross sections are displayed, calcu- 
lated for incident neutrons at laboratory energy 14.4 MeV (left panel) and 22.1 MeV (right 
panel). One may notice that a rather good agreement is obtained also in this case. Only 
at the minimum region, for 14.4 MeV neutrons, the theoretical results underestimate the 
experimental values. These small discrepancies are due to the simplicity of the MT I-III 
potential. It has been shown recently [7] that the realistic interactions further improve the 
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FIG. 2: Comparison of the momentum space AGS results (solid curves) and the current paper ones 
(dashed-dotted curves). They contain the deuteron- 12 C scattering at 30 MeV deuteron laboratory 
energy (left panel) and the proton- 13 C elastic scattering at 30.6 MeV proton laboratory energy 
(right panel). Differential cross sections for the elastic scattering and neutron stripping reaction 
are presented as a ratio to the Rutherford cross section da^/dO,. The proton analyzing power is 
displayed in the bottom of the right panel figure. The experimental data are taken from Refs. [20I — 
22\. 



TABLE II: Neutron-triton elastic (cr e ), inelastic (crj,) and total (at) scattering cross sections (in 
mb) for the selected neutron laboratory energies (in MeV) compared with the experimental data. 

E lab MT I-III Exp. 

(MeV) a e a h a t a t [Ret] 

14.4 922 11 933 978±70 [25] 

18.0 690 25 715 750±40 [25] 

22.1 512 38 550 620±24 [26] 

results, providing almost perfect agreement with the data, both for the differential as well 
as for the integrated elastic cross sections. 
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FIG. 3: Calculated n—^H elastic differential cross sections for neutrons of laboratory energy 14.4 
MeV (left panel) and 22.1 MeV (right panel) compared with experimental results of Frenje et 



ai. tea, 

Debretin et al. 



231 ] and Seagrave et al. [241 ]. 



IV. CONCLUSIONS 

A method based on the complex scaling is presented to solve the scattering problem 
above the break-up thresholds. This method does not require the implementation of highly 
non trivial boundary conditions and allows to solve the few-body scattering problem using 
square-integrable functions. Scattering problem might thus be solved using configuration 
space bound state techniques in complex arithmetics. 

Reliable and accurate results have already been obtained for the 3- and 4-body elastic 
and break-up processes, including optical potentials as well as long-range Coulomb repulsion. 
This method opens a way to explore the many-body scattering and breakup reactions and 
allows an accurate treatment of a rich variety of problems in molecular as well as in nuclear 
physics. 
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